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Entanglement of any pure state of an N x N bi-partite quantum system may be characterized 
by the vector of coefficients arising by its Schmidt decomposition. We analyze various measures 
of entanglement derived from the generalized entropies of the vector of Schmidt coefficients. For 
N > 3 they generate different ordering in the set of pure states and for some states their ordering 
depends on the measure of entanglement used. This odd-looking property is acceptable, since these 
incomparable states cannot be transformed to each other with unit efficiency by any local operation. 
| In analogy to special relativity the set of pure states equivalent under local unitaries has a causal 

structure so that at each point the set splits into three parts: the 'Future', the 'Past' and the set of 
noncomparable states. 
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I. INTRODUCTION 



Due to growing interest in quantum information the properties of quantum entanglement became a subject of 
intensive research. Entanglement gives rise to the peculiarly quantum mechanical correlations that may exist between 
| two subsystems. In spite of several papers published yearly on this subject we do not know the answers to even some 
very basic questions in this field. An example of such a question: 'Is a given quantum state (pure or mixed) of a 
composite system separable or entangled?' For the simplest case of 2 x 2 and 3x2 systems this problem is solved 
by the Peres-Horodeccy criterion which offers a simple condition sufficient and necessary for separability. On 

the other hand, in the general case of a K X M bi-partite system such a constructive condition remains still unknown 

o B8- 

Another important issue is the question: 'How much is a given quantum state entangled?'. In other words we are 
\ looking for a function of the states that quantifies entanglement. Since we do not quite understand what entanglement 
means this is a difficult question, but what we can do is to write down conditions that such a measure has to satisfy. 
A key condition is that the expected entanglement cannot increase unless the two subsystems interact in a quantum 
mechanical fashion. We therefore require monotonicity under local operations and classical communication (LOCC) 
H ||. The idea is that we have two subsystems A and B and two experimenters Alice and Bob who are allowed 
to manipulate one subsystem each in any way they please. They are also allowed to tell each other, using classical 
communication, what unitary transformations and measurements they perform on their respective subsystems so that 
Alice can allow her manipulations to depend on what Bob does and what the results of his measurements are (and 
conversely for Bob). There exists a very beautiful characterization of LOCC operations due to Nielsen which we 
will quote in due course. 

Now it is known Q that monotonicity under LOCC operations does not suffice to single out a unique measure of 
entanglement. A major point of the present paper is precisely to study this non-uniqueness in detail. The situation 
changes if we consider the asymptotic regime, in which it is assumed that we have a source producing infinitely many 
copies of our quantum system. In this situation suitable continuity requirements can be added to the list of properties 
that an entanglement measure should satisfy, and there are uniqueness theorems || that show that on pure states the 
only measure that fullfills all the axioms must agree with a particular measure known as the entropy of entanglement 
when evaluated on pure states. In the present paper we stay in the regime of one copy however. 

Classifying measures of entanglement and selecting the physically relevant ones is a subject of vivid recent interest 
[p|,p| JT^JT5[ | . One line of research starts with the investigation of invariants of local unitary transformations, like 
Schmidt coefficients for pure states |]l6| . Several such invariants for mixed states of a bi-partite system were recently 
found and 

one may try to construct existing measures of entanglement out of them or define new measures 
with attractive properties. Another possible approach is to quantify the entanglement by the distance of the given 
state to the closest separable state [^0|p^] , or to look for the best separable approximation of an entangled state p5) . 
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The aim of this paper is to elucidate some geometrical properties of the entanglement of pure states of a composite 
N x N system. Using the concept of the Schmidt decomposition we illustrate the observation ||] that in the regime 
of one copy there exists infinitely many measures of entanglement which generate different orderings of pure states. 
We establish a link between different measures based on the generalized entropies of the vector of Schmidt coefficients 
and the various distances between closest separable states 110] . 

Applying Nielsen's characterization of LOCC, that is using the concept of majorization j?],^6| (to be introduced 
below), we show that any state \tp) gives rise to a natural decomposition of the set of all pure states into four sets: 
a) the measure zero set of interconvertible states which may be obtained from by local unitary transformation, b) 
the states accessible from \ifj) by nonunitary local transformations, c) the states which can be transformed locally into 
and d) — the set of incomparable states. Once all pure states that can be transformed into each other by local 
unitary transformations have been identified with each other, the decomposition consists of only three sets, b), c) and 
d) . This picture mimics the well known structure of the light cone in special relativity, which divides the space-time 
into 'Future', 'Past' and the set of space-like, incomparable events. Hence our title (see also p7|). The same analogy 
works for the entirely different problem of non-unitary evolution of mixed states under the action of external random 
fields. Interestingly, in this case the direction of the arrow of time is reversed. 

The paper is organized as follows. In section II we review the definitions of the Schmidt decomposition of a pure 
state, quantum Renyi entropies and majorization, while in Section III we discuss different measures of entanglement, 
analyze nonhermitian evolution of the density matrices, for which the majorization theory is applicable, and the 
corresponding local transformations of pure states. Propositions emphasizing the geometrical interpretation of the 
Weyl chamber consisting of ordered eigenvalues of a density matrix (or ordered vector of the Schmidt coefficients for 
a pure state of a bi-partite system) are proved in Appendix A. 

II. MATHEMATICAL TOOLS AND DEFINITIONS 

A. Schmidt decomposition 

Consider a pure state of a composite Hilbert space TL — TLa^'Hb of size N 2 . Introducing an orthonormal basis 
in each subsystem, we may represent the state as 

N N 

n) ® \m). (2.1) 

n— 1 m— 1 

The complex matrix of coefficients C of size N needs not to be Hermitian nor normal. Its singular values (i.e. square 
roots of eigenvalues \k of the positive matrix C^C) determine the Schmidt decomposition f28| , |l6| , p9f 

N 

W) =X!\/A#/><g>|fc"), (2.2) 
fe=i 

where the basis in Tt is transformed by a local unitary transformation U <8> V. Thus \ki) = U\k), and \kn) — V\k), 
where U and V are the matrices of eigenvectors of C^C and CC\ respectively. The remarkable thing is that only a 
single sum is involved. In the generic case of a non-degenerate vector A, the Schmidt decomposition is unique up to 
two unitary diagonal matrices, up to which the matrices of eigenvectors U and V are determined. The normalization 
condition = 1 enforces X)fe=i = 1- Thus the vector A = (Ai, Xn) lives in the (N — 1) dimensional simplex 

Spf. Note that the Schmidt coefficients A*, do not depend on the initial basis |n) ® |m), in which the analyzed state 
is represented. They are invariant with respect to any iocai operations U a ® Ub, and thus they may serve as 
ingredients of any measure of entanglement. The Schmidt simplex is almost but not quite the same as the space 
of orbits under local unitary transformations; since we can change the ordering of the eigenvalues by means of local 
unitaries some further identification of points in the Schmidt simplex has to be done before we have the space of 
orbits. This is discussed below (section HID). 

The Schmidt coefficients of a pure state \ip) are equal to the eigenvalues of the reduced density operator, obtained 
by partial tracing, p A — tr b ■ The pure state is called entangled if it can not be represented in the product 
form \ip) = \iPa) ® \iPb), where \ipA) € Ha and € 7~Cb- This is the case if and only if there exists only one 

non-zero Schmidt coefficient, Ai = 1. 
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B. Entangled mixed states 



Mixed states p = Pi\(fi) (<fi\ with pi > and YliPi = 1 will also be useful in our considerations. A state p is 
called a product state if it can be represented as a tensor product, p pro d — p A <£> p B , where p A acts in Ha and p B 
acts in TIb- A mixed state p is called separable if it can be represented as a convex combination of product states, 
Psep = J2j IjPf ® Pj ') where > and ^ ■ = 1. A mixed state which is not separable is called entangled. It is 
easy to see that for pure states both definitions are consistent. 



C. Measures of entanglement 

There exist several different possibilities to quantify quantum entanglement - for a review on this subject see e.g. 
J|,[l5]]. Following Vedral and Plenio |nj we assume that each measure of entanglement E(p) has to satisfy the following 
conditions: 

(El) E(p) = if p is separable, (the condition 'if and only if occurs to be too strong |l5[| ) 
(E2) Local unitary operations do not change the entanglement, i.e. E(p) = E(Ua <8> UspU A (8 U B ) 
(E3) The expected entanglement cannot increase under operations involving local measurements and classical 
communication (LOCC), followed by post-selection: 

M 

J2tr Pi E(-*-) <E(p), (2.3) 
i ^ trpi / 

where pi = VipV} , each of the operators Vi is local, i.e. Vi — V { A ® V^ , and the set of M these operators defines 
an positive operator valued measure (POVM) H|, i.e. Y^tLi v i v i = 1 Note that the definition allows the actions 
on the two subsystems to be correlated; this is how classical communication enters. Note also that it is the expected 
entanglement that cannot increase — if we end up with a statistical ensemble of states there may be a non-zero 
probability of increased entanglement. 

Sometimes one imposes further requirements 

(E4) Additivity: E( Pl ® p 2 ) = E(p x ) + E(p 2 ), 

(E5) Continuity: E is a continuous function of /?, 

or 

(E5') Asymptotic continuity: E is a continuous function of the fidelity for n copies of the same pure state, \ip)® in 
the asymptotic limit n — > oo, 

the necessity of which is still disputed in the literature |8|,|l5|,H ■ Condition (E4) is most welcome but difficult to prove 
in the general case of arbitrary density matrices, while condition (E5) is satisfied by many different measures. On the 
other hand, condition (E5') is rather strong: assuming that the measure of entanglement fulfills some technical condi- 
tions which quantify asymptotic continuity one may show |||n| that for pure states it is proportional to the entropy 
of entanglement, E(\ij)){ip\) = iJjv(IV')) = — tr(p^) ln(p A ). Another set of axioms for the measure of entanglement 
which singles out the entropy of entanglement is given by Rudolph Such requirements are indeed appropriate 

in the asymptotic regime but it remains interesting to investigate the regime of one copy where only the first three 
conditions are imposed. One of the aims of this paper is to emphasize that there exist several reasonable measures of 
entanglement, which for pure states arc different from the entropy of entanglement. 



D. Quantum Renyi entropies 



Consider an N dimensional vector x with non-negative components normalized as Yli=i x% = 1. The distribution 
of Xi may be described by the Shannon (information) entropy Hn '■= — ^2^ =1 Xilnxi. A more general family of 
quantities characterizing the components of x is provided by the Renyi entropies [j3l[ often used in information theory 

©; 

i N 

H a (x) -j— 1*[5>"]- (2-4) 

i=l 
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As usual, in the definition of entropies we adopt the convention that OlnO = 0, if necessary. The non-negative number 
a =/= 1 is a free parameter labeling the generalized entropy. It is easy to see that lim a ^i H a = Hn, so for consistency 
we will write sometimes Hi for H^. The quantities H a (x) vary from zero (for x with one non-zero component) to 
In TV, for the vector with all components equal, Xi — 1/TV. It is possible to show that the generalized entropy is a 
non-increasing function of its parameter: for any x and a.2 > cti, the inequality H a2 (x) < H ai (x) holds |3^] . 

Some special cases of H a are of particular interest. Let r denote the number of non-zero components of the vector 
x, x\ be its largest component, and \x\ = (X)i=i x l be its length. Then 

H (x)=lnr, H 2 {x) = -\n\x\ 2 , H^x) = -\nx x , (2.5) 




(001) (001) 




(100) (010) (100) (010) 

FIG. 1. Iso-entropy curves H a (x) in the space of TV = 3 vectors x plotted for (a) a = 1/3, (b) a = 1 (also contour lines of 
the Shannon entropy), (c) a = 2 (circles of points equidistant from a?*), and (d) a = 5. 



To get some experience with the properties of generalized entropies let us look at Fig. |l|, which shows the iso-entropy 
curves for TV = 3. The equilateral triangle, of side length y/2, lies in the plane x$ — 1 — xi — x\ and its corners are 
labeled by the co-ordinates (x%, x 2 , X3). At these three points the generalized entropies attain their minima, H a = 0, 
for any values of a > 0. The maximum H a — In 3 is achieved at the center of the triangle, which represents the 
uniform vector x*. The maximum is rather flat for a — 1/3, as shown in Fig. |l|a. For even smaller, non-negative 
values of a the isoentropy lines become parallel to the closest side of the triangle. The entropy Ho vanish at the 
corners of the triangle, is equal to In 2 at its sides and equals to In 3 for any point inside the triangle. The other 
example, a = 5, presented in Fig. [l]d. is similar to the limiting case H^, for which the iso-entropy curves are parallel 
to the most distant side of the triangle. 

For any two values of a the isoentropy lines intersect. Thus any pair of two different generalized entropies induces 
a different ordering of vectors x. In our considerations the vector x will represent two entirely different objects: 

a) the spectrum d of a TV = 3 density matrix p 
or 

b) the vector of the Schmidt coefficients A of a pure state of the 3x3 composite system. 

Thus the pictures of isoentropy lines may have two entirely different meanings. In the former case a), the entropy 
-Hat is equivalent to the von Neumann entropy, Sn(p) = ~Tr[plnp], while other quantum Renyi entropies 

S a (p) = — !— lntrp (2.6) 
1 — a 

play the role of different measures of degree of mixing (3^]34|] . The entropy Sq is a function of the rank r of the density 
matrix, i.e. the number of positive eigenvalues of p. Thus it is not a continuous function of its argument in contrast to 
generalized entropies S a with a > 0. For a — 2 we have S2 = — ln[Trp 2 ] = In R, where the purity (also called inverse 
participation ratio) R(p) — [Trp 2 ]" 1 describes the "effective number of states" involved in the mixture p and varies 
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from unity (for pure states) to N (for =1/3). The purity of any state depends only on its Hilbert-Schmidt distance 
from the maximally mixed state p„, represented by the center of the triangle. Thus the curves of constant R have the 
shape of a circle, and the corners - the points most dis tan t form p» - denote three orthogonal pure states. Note that 



due to the presence of the logarithm in the definition (2.6) the generalized entropies are additive for product states, 



S a (pi ® pi) = S a {pi) + S a (p 2 ) (2.7) 

for any a > 0. 

In the latter case b) the centers of the triangles plotted in Fig. 1 represent the maximally entangled states |i/>*} = 
(|11) + |22> + |33))/V3 and the corners denote separable states. The generalized entropy of a vector of the Schmidt 
coefficients equals the Renyi entropy of the reduced density matrices, H a (\ip)) — S a [trB{\ip) — *S l a[t I 'A(|'0)(V ; l)]- 
In particular Hn coincides with the entropy of entanglement of the corresponding pure states of the composite 3x3 
system. As discussed later, all generalized entropies H a for a > fulfill the conditions (E1)-(E3), and may thus serve 
as legitimate measures of the pure state entanglement. Moreover, for any prod uct pure states these measures are 
additive, H a (\ipi) ® IV^)) = H a (\ipi}) + H a {\ip2))- This property following from fl2.7| ), was pointed out by Vidal ||. 

Note that an analogous construction of isoentropy hypersurfaces in the (N — 1) dimensional simplex may represent 
mixed states acting in Hn or pure states of the composite N x N system. For any N > 3 the entropies H a (x) provide 
not equivalent measures of uniformity/disorder in the distribution of the components x%. 



E. Majorization 



The question, which of any two vectors is more 'mixed' (or more disordered) can be approached by the theory 
of majorization |35|, which allows one to introduce a partial order into the this set. Consider two vectors x and y, 
consisting of N non-negative components each. We order the components of both vectors in a decreasing order, (what is 

sometimes emphasized by the notation x^), and assume that they are normalized in the sense Xi = 5~)j—x Vi = 1- 

We say that x is majorized by y, written x -< y , if 



2 Xi - £ Vi 



(2.8) 



for fc=l,2,...,JV— 1. Since the sum of all components of both vectors are equal, we obtain an equivalent formulation 
of majorization 



JV N 
i=l 



(2.9) 



i=l 



for all I = 2, . . . , JV. Vaguely speaking, the vector x is more 'mixed' than the vector, y and its distribution function 
grows slower with the index i (see Fig. 0). 




1 2 3 4 5 6 7 



FIG. 2. Idea of majorization: (a) the ordered vector = {xi, . . . , £10} (•) is majorized by y = {yi, . . . , yio} (A), since the 
corresponding distribution function takes larger values (b), and both functions do not cross. 
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The functions which preserve the majorization order are called Schur convex, 



x<y implies f(x) < f(y). (2.10) 

Examples of Schur convex functions include f q (x) — x t f° r an y <7 > 1 , fq{%) — ~ X)i=i x 1 f° r an y < g < 1, and 

} e {x) — ^2iL± x i hi a;, j35j, while these functions with reversed sign are Schur concave. To show a simple application 
of the theory of majorization let us then consider two density matrices p x and p y , with spectra x and y, respectively. 



If x -< y then due to Schur-convexity H a (x) > H a (y) for a > 0, since at a = 1 the prefactor 1/(1 — a) in (2.4) changes 
the sign and the direction of the inequality is reversed. 

A classical theorem by Hardy, Littlewood and Polya states that x -< y if and only if there exists a bistochastic 
matrix B such that x — By. A matrix By with non-negative elements is called bistochastic if sum of all elements 
in each row and each column equals to unity. As shown later by Horn p6[ this theorem can be strengthened by 
replacing the word bistochastic by unistochastic. A bistochastic matrix B is called unistochastic, if there exists an 
unitary matrix U such that By = \Uij\ 2 . All N — 2 bistochastic matrices are unistochastic, but this is not true for 
7V>3|||. 

Another useful feature of majorization is related to the so called T-transform. It is a matrix of size N which acts 
as identity on all but two dimensions where it has the form of a bistochastic matrix 



T 



t l-t 
1-t t 



(2.11) 



so t G [0, 1]. One can prove p5| , p7[ that x -< y if and only if there exist a set of T-transforms Ti,T%, ...,Tk that 
x = T1T2 ■ ■ ■ Tky. The maximal number of the T-transforms needed is not larger than N — 1. A simple proof of this 
fact and of the Horn lemma is recently provided by Nielsen [M. 

Sometimes it is useful to compare non normalized vectors, such that 5Z i=1 Xi ^ X)i=i Vi ~ 1- I n such a case one 



introduces the relation of weak majorization |35| defining x -< w y if 



k k 

^Xi<^ yi for k = l,2,...,N. (2.12) 

i=l i=l 

Then vector x is said to be weakly submajorized by y, or alternatively, weakly majorized from below. 

Majorization technique allows one to obtain effective necessary criteria for separability of mixed states, as recently 
shown by Nielsen and Kempe |39f| : If a mixed state p is separable then d -< cIa and d^ds where dA and ds denote 
the spectra of the reduced operators p A and p B . Since the Renyi entropies are Schur concave, one obtains 

S a (p A ) < S a (p) and S a (p B ) < S a (p) (2.13) 

for any separable p and a > 0, which generalizes earlier results of Horodeccy Q. 



III. PURE STATES ENTANGLEMENT 

A. Distance to the closest separable state 

For any entangled state p en t the distance to the closest separable state p c may be considered as a measure of 
entanglement. Each metric in the space of mixed quantum states will thus define a certain measure of entanglement. 
Vedral and Plenio have shown (To) that the Bures distance Db(pi, P2) 2 = 2 — 2tr| y o^ 2 y 02y0 1 ^ 2 1 , to the closest 
sep arab le state fulfills the required conditions (El) — (E3). For a pure state written in the Schmidt decomposition as 



in ( 2.2) the closest separable state is 



N 



= J2^k\kk){kk\ , (3.1) 



fe=i 



(this result is proved in |T(| for N = 2 and stated for larger N) . The state p c is typically mixed and its squared Bures 
distance to the pure state p^ = reads 



£>s(PV»Pc) = 2-2 



\ 



^A2=2-2v^7i?. (3.2) 



fc=i 



G 



The inverse participation ratio R is closely related with the generalized entropy of order 2, namely R = exp(i?2)- 
Therefore the circles around the maximally entangled state, forming the isoentropy lines for H2 (see Fig. lc), form 
also the lines of the same Bures entanglement. The same concerns generalized concurrence, another measure of 
entanglement recently introduced by Rungta et al. [ (tl| . 

Vedral and Plenio discussed also the quantum relative entropy, S(p±\p2) '■= tr[/9i(lnpi — lnpa)]- Although this 
quantity does not induce a true metric (e.g. it is not symmetric), the entropy S(p en t\ p c ) satisfies (El) — (E3) and may 
be considered as a measure of entanglement. Moreover, they showed that the state ( |3.l[ ) is also the 'closest' separable 
state to the pure state p^p and the smallest relative entropy SWun = J3fe=i ^ k m ^ k coincides with the Shannon entropy 
Hn of the vector of the Schmidt coefficients. This quantity, often briefly called entanglement of the pure state p^,, 
is given a special physical meaning, since the probability of not distinguishing between the analyzed entangled state 
p en t and the closest separable state p c after n measurements behaves as exp[— nS(p en t\Pc)] [P-0| - Furthermore, the 
entanglement of formation of any mixed state p, defined as the minimal average entanglement of pure states, 

the mixture of which generates p, for a pure state reduces to the entropy of entanglement Hn . 

For any entangled pure state (2.2) one may also look for the closest separable pure state. This problem was 
recently considered by Lockhart and Steiner Ji^| , who show that the state \<p) := corresponding to 

the largest Schmidt coefficient Ai :— max{A;} is the closest. They were using the Hilbert-Schmidt distance, 
Dhs(pi, P2) — [tr(pi — P2) 2 ] 1 / 2 , and found DnsiPip, P&) = [2 — 2A1] 1 / 2 . However, the projective cross-ratio 
k := \(ip\(j>}\ 2 — Ai is the only parameter determining the standard distances in the space of pure state, e.g. the 
Bures distance Z3s(p^,p^) = [2(1 — ^JTl)] 1 / 2 , the trace distance D tr (p^,p<p) — 2[1 — k] 1 / 2 , and the Fubini-Study 
distance Dpsd^), |</>)) = iarccos(2K — 1) = arccos( v / K). All these functions are monotone, so the measure of entan- 
glement defined as the distance (any of the above) to the closest separable pure state defines an ordering of the pure 
states identical with that given by the Chebyshev-like like entropy = — lnAi. We have thus shown that three 
different settings of the problem of finding for any pure state of an N x N system the closest separable state reduce 
to the generalized entropies with a = 1, 2 and 00. 

On the other hand one should not expect that every reasonable mea sure of pure states entanglement satisfying 
conditions (E1)-(E3) is neccessarily a function of one of the entropies (EM). As a counterexample let us mention 
the coefficients of the characteristic polynomial of the nontrivial block of the Gram matrix introduced in |2^] . They 
can be expressed as symmetric functions of all Schmidt coefficients, e.g. J2i>j^i^j' Y^,i>j>i A-iAjA;,..., where the 
summation goes over all possible sets of k indices, k — 2, ...,N [Q. All these functions are Schur-concave |5f|, and 
thus entanglement monotones, although in the general case, (for N > 3) they are not functions of the generalized 
entropies. 



B. 2 x 2 system - only one ordering of pure states 

For pedagogical reasons we start to analyze the pure state entanglement with the simplest 2x2 system. Pure states 
of this N — 4 system may be parametrized as 

\ip 4 ) = (cos^^in^cos^e^^in^sin^cos^ie^^in^sin^sin^ie^ 1 ), (3-3) 

where S [0, 7r/2], and ipu G [0, 2ir) for k = 1,2,3. The states belong to the complex projective manifold CP 3 . 
This space is compact and has 6 real dimensions, e.g. three polar angles $i and three azimuthal angles tfi. 

To get a better understanding of CP 3 , let us recall the structure of the space of pure states for systems of even lower 
dimension. The set of the N = 2 pure states, ^2) = (cost?, sin$e^), is described by the Bloch sphere CP 1 ~ S 2 . 
The sphere may be drawn in a simplified way by an interval (a meridian labeled by $ G [0,7r]), each point of which 
represents a circle (a parallel 1? =const with <p S [0, 2ir)). At both poles the circles reduce to a point - see Fig. 2. 

The 4-dimensional manifold CP 2 of the N — 3 pure states can be parametrized as | ^3 ) = 
(cos^^in^costJie^^in^sin^ie^ 1 ) where $i S [0,tt/2], tf 2 6 [0,tt/2) and (p 1 ,tp 2 £ [0,2tt). These local coor- 
dinates allow us to describe (almost all of) this space as shown in Fig. 2b. The angles ($1, # 2 ), describe a point in the 
positive octant of a sphere 5 12 (or in an equilateral triangle, what is topologically equivalent) which represents entire 
2-torus formed of both phases ipi jlBj . Each point on one the three edges of the octant represents a circle, so each 
entire edge corresponds to a sphere. Note that three corners of the triangle are not the 'corners' in CP 2 - in the same 
sense as for N — 2 the poles (the edges of the meridian) are topologically equivalent to all other points on the sphere. 
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a) N = 2 b) N = 3 




s cp'=s 2 CP 2 



FIG. 3. The Bloch sphere, S 2 = CP 1 , of all pure states of a two levels system (a) may be drawn schematically as a line, 
each point of which represents a circle (c). In the same, simplified manner we depict CP 2 - the set of all pure states of the 
N = 3 system (b). Each point inside the octant of a two-sphere, associated with the angles ($i,i?2), represents a torus T 2 
spanned by the phases (951, ipi)- Each point at the edges of the octant denotes a circle, so each of three edges represents a 
sphere. The questionmark represents the 4 dimensional manifold CP 2 that we are not able to reproduce exactly in the picture. 



These two examples help us to imagine the structure of the space CP 3 of the N = 4 pure states. It may thus be 
represented as a 1/16 part of the hypersphere S 3 , parametrized by the angles Each point inside such a 'hyper- 
octant' (or simpler, the tetrahedron) represents a 3-torus spanned by the phases tp^ each point on the face of the 
simplex a 2-torus and each point of the edge a circle. In general all points of CP 3 are equivalent. This symmetry 
becomes broken if out of these N = 4 system we distinguish two two-levels subsystems. 

The corners of the tetrahedron represent mutually orthogonal states; let us use the standard notation and label 

them ( ), ( — h), (H — ), and (H — h). In principle one could find analytically the entropy of entanglement H\ (or any 

other generalized entropy H a ) for each point of the simplex (and each choice of the phases tfk)- To get a more 
transparent picture we prefer to plot Hi(\ip)) at the surface of the tetrahedron - see Fig. 4. 

Four corners represent the separable states. The same concerns four edges of the simplex irrespective of the phases 

running along the circles. Find two maximally entangled Bell states proportional to | H — ) + | — h) and | + +) + | ) 

localized in the middle of two non-connected edges. The other two orthogonal entangled states, (| H — ) + | — h))/\/2 

and (| + +) — I ))/v2, contain non-zero phase, (minus sign) and do not belong to this plot. The structure of the 

picture changes with other choices of the basis. If one uses the Bell basis of the four maximally entangled states and 
places them into four corners then the separability of the edge states depends on the phases (p. 

One may portray any other generalized entropy H a (\ip)) in a similar plot. However, in this case of pure states of 
the 2x2 system, all measures of the entanglement are equivalent in the sense that one measure is a function of the 



other one. Therefore any two entanglement measures, E\ and E 2 , generate the same ordering 46 



E 1 {p 1 ) > E 1 {p 2 ) o E 2 {p x ) > E 2 (p 2 ) (3.4) 

for any density operators representing pure states, pi = \ipi)(ipi\. This is due to the fact that for N = 2 the 
entanglement of any pure state is completely characterized by only one relevant parameter - the Schmidt angle 
(3 € [0, 7r/4], such that A1 = cos 2 (/3) while X 2 — sin 2 (/3). The distribution of the Schmidt angle for random pure 
states distributed according to the natural, unitarily invariant measure on CP 3 is given by P(/3) = 3cos(2/3) sin(4/3). 
A simple integration allows us to find the mean angle, (/3)cp 3 = 1/3 which incidentally equals the mean entropy of 
entanglement, (iTi(|^))) CP 3 = 1/3 @>@. 
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(— ) 




( — ) (++) C— ) 



FIG. 4. Entropy of the entanglement, H\{\ip)) for pure states of the 2x2 problem at the surface of the t?-tetrahedron 
representing CP 3 (dark colors denote large entanglement). Labels refer to the corners representing the standard basis. Copy 
this picture magnifying, cut out the triangle, fold three times along the edges and glue together into a thetrahedron to enjoy 
the symmetry of this representation. 



C. Global transformations of mixed states 



Before dealing with the more complicated case of entanglement for the 3x3 system let us make a detour to have 
a closer look at the space of mixed quantum states acting in TLn- Any state p may brought to a diagonal form by a 
unitary evolution p — > p' = UpU* which preserves the vector of eigenvalues d = {di,d2, ...,djv}- The normalization 
condition assures trp — Y2i=i di = l. 

Any physical operation can be described by a superoperator M : p — > p' that preserves positivity of the mixed 
state p. Strictly speaking the term 'positive' refers to positive semidehnite Hermitian operators, which do not have 
negative eigenvalues. Moreover, the system may be coupled to an environment, and M may be trivially extended to 
M <g> I. A superoperator M is completely positive if for all such extensions M ® I is positive. A completely positive 
(CP)-map which preserves the trace (normalization of the state is conserved) is called a stochastic map or quantum 
channel. It is capable to describe any quantum operation including the process of quantum measurement and may 
be represented in the Krauss form Q 

k 

p^ p' = M{p)=Y J A lP A\, (3.5) 

i=l 

where the set of k operators Aj satisfies the completeness relation, ^i=i A]Ai = I. 

In the following we shall consider a smaller subset of quantum channels, which may be written in the form p' = 
M(p) = ^2i=iPiUipUj , where each operator Ui is unitary and the positive coefficients pi sum to unity. These 
operations, called external random fields ]5Q], are unital, since the maximally mixed state p* = I/N is preserved, 
M(p*) = p*. Maps which preserve both the trace and the identity are called bistochastic. 

Let d' denote the vector of eigenvalues of p' = M (p). Then one can prove J5l],|5^] the following majorization relation, 
a -< d. Using the concept of the T-transforms we are going to show the converse: given a vector of eigenvalues d! such 
that d' -n d there exists a bistochastic map M : p' — M(p), written p gl 2 1 ^ 1 p 1 Thus both properties are equivalent, 







(3.6) 



For simplicity we start discussing the N = 3 case, represented in Fig. 5. Consider a state p with spectrum d, (its 
components are ordered as d\ > > cfe) and a transformation M%2(p) — wnUnpU^ + (1 — w\2)p. The unitary 
matrix [/ 12 acts as an identity on all but two first components, (labeling the matrix) which get mixed by the orthogonal 
submatrix: (JJ 12)11 = (^12)12 = (£^12)22 = l/V^ and (JJ 12)21 = — 1/V2 ■ In this way Mi moves p along the horizontal 
line joining d with d a = {d\2, di2, cfe), where d\2 — (d% + d2)/2. The length of this move is controlled by the weight 
parameter w\2- An analogous transformation Af23(p) defined by the unitary matrix U23 which preserves the first 
eigenvalue d\, allows one to obtain any state along the line between d and db = (d\, ^23, ^23), where ^3 = (c?2 +da)/2. 
It is easy to see that any state with the spectrum fulfilling d! -< d (shaded region in Fig. 5a), may be reached by 
a superposition of operations Mij, which correspond to T-transforms. Observe that the lines limiting the accessible 
region are parallel to the isoentropy curves Hq and i?oo, respectively. 




(100) {V2V2Q) (100) (010) (100) (010) 



FIG. 5. Simplex of eigenvalues for N = 3 density matrices: corners represent pure states, center the maximally mixed state 
p* = 1/3; (a) asymmetric part of the simplex (Weyl chamber), for which di > di > dg. Evolution governed by two T-transforms 
defines the shaded region accessible from p; (b) the shape of the 'light cone' depends on the degeneracy of the spectrum: F - 
denotes Future, P Past, and C the non comparable states; (c) splitting of the entire simplex into three parts with permutation 
operation allowed. 



This construction is simply generalized for arbitrary TV: the accessible region of states with spectra majorized by d is 
a convex polytope with 2(N — 1) faces. The point d sits in one of its corners, which joins N — 1 edges - lines generated 
from p by the transformations M^+i, where i = 1, ...,7V — 1. Therefore every face forming the part of boundary 
of the accessible region, (corresponding to the light cone in special relativity), may be generated by a sequence of 
(N — 2) T-transforms acting on d. The first transform, M\ t 2, moves p along a line paralell to the closest face of the 
simplex, so the entropy Hq is conserved. All other transforms, M^+i, do not influence the largest component, d\, so 
the entropy remains unchanged. 

The majorization property d' -< d assures that the Von Neumann entropy and all other entropies S a for a > 
grow (strictly speaking do not decreas e) d uring the time evolution governed by random external fields M. Note that 
the factor 1/(1 — a) in the definition (2^4) is crucial to assure the Schur convexity of H a . For any state p it is easy 
to specify the set of all states p' accessible from it by external random fields. For these states, denoted in Fig. 5b 
by F as 'Future', all generalized entropies satisfy S a (p') > S a (p). The complementary set of all states 1 p which may 
be transformed by M into p is denoted by P as 'Past'. Obviously S a (p) > S a ('p). There exists also a set C of 
noncomparable states, which are not ordered with respect to p by the majorization relation. This structure is thus 
analogous to the special relativity picture of the light cone. The actual shape of the 'cone' (e.g. the size of the angle 
'Future') depends on the degeneracy of d as shown in Fig. 5b: in the generic case it is equal to 27r/3, while at a 
bisetrix it shrinks to 7r/3. The 'world line' crossing the point p and representing a possible trajectory in the space of 
density matrices is located inside the accessible region F. 

The simplex can be divided into N\ asymmetric simplices (called Weyl chambers in group theory), which differ by 
the permutation of the components of the vector d. In the discussed case N = 3 three bisetrices split the simplex into 
6 right triangles, one of which is shown in Fig. 5a. A permutation of the components of d, which might be considered 
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as a special unitary case of the random fields M, maps d into a corresponding symmetric point in one of the other 
6 right triangles. Therefore, in the generic case of a non-degenerate spectrum of p, the entire decomposition of the 
N = 3 equilateral triangle of eigenvalues looks as shown in Fig 5c. The set 'Future' has the form of a hexagon, which 
contains two triples of equal sides, and becomes regular for d = (3, 2, l)/6, or in general for d = (1/3 + x, 1/3, 1/3 — x), 
where x € (0, 1/3]. If the spectrum of p is degenerate (the vector d is located at a bisetrix), the hexagon reduces to 
an equilateral triangle. 

It is not difficult to understand how this construction works in higher dimensions: the set 'Future' is a convex 
polytope, defined by the convex hull of N\ points obtained from d by permutations of its components. This set 
contains N\ identical simplices surrounding and is topologically equivalent to a ball. It consists of Nl vertices (one 
for each permutation), each joining N — 1 edges. Thus in the generic case this polytope has N\(N — l)/2 edges and 
2^ — 2 (hyper)faces of dimension N — 2. For N = 4 it is a polyhedron with 24 vertices, 36 edges and 14 faces; if the 
vertex at d lies on the line 

d(» = |(3,2, 1,0) + i^£(i, 1,1,1); 0<x<l (3.7) 

then we obtain an Archimedean solid known as the truncated octahedron whose surface is composed of 8 regular 
hexagons and 6 squares. By definition an Archimedean solid has regular but not equal faces |^3| . If the eigenvalues 
are degenerate, the polytope degenerates as well. It is amusing to notice that a number of Platonic and Archimedean 
solids appear for special choices of d. In particular there is a surface in the simplex at which the degeneracy pattern 
is d = (x , x, x, X3). If the vertex sits at a particular line on this surface then the polytope is the Archimedean 
cub-octahedron. Similarly, the surfaces determined by d = (xq, x±, x, x) or d — [x, x, X2,%a) contain lines giving rise 
to the truncated tetrahedron. On the lines d — (xo, x, x, x) and d = (x, x, x, X3) we obtain regular tetrahedra while the 
line d = (x,x,y,y) yields regular octahedra. So much about the 'Future'. The set 'Past' consists of TV disjoint parts 
occupying neighbourhoods of each corner, while C consists of 2 N — 2 parts, each for every face of the set 'Future'. 
The flat geometry of the simplex of eigenvalues has a concrete meaning, since the Euclidean distance between any two 
points of a Weyl chamber (Fig. 5a) is equal to the Hilbert-Schmidt distance between the global orbits they generate 
- see Appendix A. 

Note that the analogy with special relativity does not work in the simplest case of N — 2, for which all measures of 
the degree of mixing are equivalent. All of them become functions of one single parameter, (say, the distance to the 
center of the Bloch ball), and generate the same order of the density matrices. The situation changes for N = 3, for 
which the spectrum of p depends on two parameters, as recently discussed by Jozsa and Schlienz ]54| . They obtained 
interesting relations between different entropies H a , which may also be proved using the techniques of majorization. 
which may also be proved using the techniques of majorization. For example, if we move in the triangle of eigenvalues 
along a circle of a constant R and trp 3 grows (the entropy H3 decreases), then the von Neumann entropy Hi increases. 
Such relations may be figured out by superimposing together different isoentropy curves shown in Fig. 1. 



D. Local transformations of pure states 

Consider now the set of pure states of an N x N composite system. We divide it into equivalence classes consisting of 
states that can be transformed into each other by local unitary transformations; the resulting orbit space is precisely a 
Weyl chamber in the Schmidt simplex. Let A^, = {\i(ip), Ajv(V')} represents the vector of the Schmidt coefficients 
of a state ordered in the descending order. Define a deterministic LOCC transformation as one that takes states 
to states (as opposed to one that transforms a state to a statistical ensemble of states). Nielsen has shown that a 
deterministic LOCC transformation of a given state into \4>) is possible if and only if the corresponding vectors of 
Schmidt coefficients satisfy the majorization relation (?]] 

\^)^\4>) A^A . (3.8) 

Note a similarity between this relation concerning the local transformation of pure states and the condition ( |3.6| ) 
concerning the global bistochastic transformations of mixed states. 

Any pure state \ip) allows one to split the Schmidt simplex, representing local orbits of pure states, into three 
regions: 'Future' F, 'Past' P, and non-comparable, C. This structure for N = 3 is shown in Fig. 6. Although Fig. 
5 describes an entirely different physical problem, both picures are almost identical. The only difference consists in 
the arrow of time: the 'Past' for the evolution in the space of density matrices corresponds to the 'Future' for the 
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local entanglement transformations and vice versa. The triangles plotted here should not be confused with Fig. 4. In 

fact the analogous plot for N — 2 would contain only one line, say joining the states (++) and ( ), along which 

the entanglement takes all values in [0, ln2]. On the other hand, the triangle shown in Fig. 6 may be given a direct 
geometrical interpretation, provided it is deformed into an octant of the sphere. Formally this may be achieved by 
the transformation {Ai, An} — > {A^, A^}. Then the length of the arc joining two points of the Weyl chamber 
(which covers 1/48 of the sphere) is proportional to the Fubini-Study distance between the local orbits they generate 
- see Appendix A. 

In general, all different measures of mixing increase during the operations considered and all the measures of 
entanglement decrease. In particular it follows that all generalized entropies H a (A) and their convex functions belong 
to the class of entanglement monotones which do not increase during local operations [|| . For example Hoo (A) equals 
minus maximal Schmidt coefficient, — Ai, a simple entanglement monotone discussed by Vidal The entropies H a 
satisfy thus the conditions (E1)-(E3) required for a measure of the entanglement. The maximal fidelity of a pure 
state, Fmax(V') = Ei=i V~Ai) 2 /N [H), a convex function of the generalized entropy, F max = exp(i?i/ 2 )/./V, is also 
an entanglement monotone due to the Schur convexity. So is another measure of entanglement called robustness, 
which is defined as the minimal amount of separable noise that has to be mixed with the analyzed state to wash out 
completely its quantum correlations [§6]j|]. For pure states the robustness equals NF max — 1 = cxp(iJ 1 / 2 ) — 1- 

Another quantity used to characterize the entanglement is called negativity |57],Q. It is based on the concept 
of partial transposition (transposition in one of the two subsystems). It is known that for any separable state p its 
partial transpose p T2 is positive |^,^|, and the negativity is expressed by the trace norm, t := \p T2 \ti — 1, i.e. by the 
sum of absolute values of the eigenvalues. Negativity is an entanglement monotone and satisfies the axioms (E1)-(E3) 
^9[. For pure states of 2 x 2 system it is equal to concurrence, used by Wootters to calculate the entanglement of 



formation [|42 13 1. For any N x N pure state written in the Schmidt form ( [2.2] ), the partially transposed matrix p T2 
has a simple block structure: it consists of N Schmidt numbers at the diagonal (which sum to unity) and N(N — l)/2 
blocks of size 2, one for each pair of different indices Eigenvalues of each block are i^/AjAj, so the negativity 

reads 



N N 



t(M) = 2 J2 V / ^Aj=(Ev / ^J "I. (3-9) 

i>j = l i—1 

which varies from for separable states to N — 1 for maximally entangled states. It is easy to see that for pure 
states the negativity is equal to robustness, so is also a function of the quantum Renyi entropy entropy of order one 
half, t = exp(ii!/ 2 ) — 1. Interestingly, the same behaviour for the pure states is characteristic to the cross norm - 
another measure of entanglement introduced by Rudolph E3] . Although negativity t does not satisfy (E4) its nonlinear 
function, H 1 / 2 = ln(t + 1), is additive in the general case of arbitrary mixed states, as all other Renyi entropies. On 
the other hand, the negativity suffers a serious drawback: it is not capable to detect so-called bound entangled states 
with positive partial transpose, which are known to exist for quantum systems with dim7i > 8 p^JHJ . 

Let us emphasize that the analogy to special relativity holds for the description of local transformations of pure 
states entanglement. This analogy was put forward explicitly in the paper by Hwang et al. J27j] and used (correctly) 
to describe a related problem of the mixed states entanglement. As we point out in this work, the same analogy is 
useful also in investigating the pure states entanglement. However, the authors of |27j discussed only the simplest 
2x2 problem, for which all entanglement measures for pure states generate the same ordering. In the general case 
this is no longer true. For example the maximal fidelity F max and the entropy of entanglement Hi generate different 
ordering of pure states for N > 2, as can be deduced in Fig. 1 for N = 3. The fact that various entanglement 



measures generate different ordering of the set of pure states (condition (3.4) is violated) , has been raised in ^0|.|l2[ . 
Vidal concludes that it is necessary to use N —1 different measures of entanglement to provide a detailed description 
of the set of pure states of N x TV system. This is a consequence of the fact that there exists exactly N — 1 independent 



Schmidt coefficients in the decomposition (2.2). A possible set of entanglement monotones is given by the sum of the 



N — k + 1 smallest Schmidt coefficients, Ek = Yli=k ^ with k = 2, ...N Note that the sum of N — 1 smallest 
coefficients is related to the already mentioned monotone, the largest coefficient, Ai = 1 — E 2 - 
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FIG. 6. Simplex of Schmidt coefficients for 3x3 pure states: corners represent separable states, center the maximally 
entangled state \tp*)- Panels (a-c) analogous to these in Fig. 5, observe the opposite direction of the arrow of time. As in Fig. 
5 the bounding lines are the level curves of Ho and H x familiar from Fig. 1. 

Pure states with the same set of the Schmidt coefficients (e.g. 6 points in fig. 6c generated be reflections with 
respect to all bisetrices) are called inter convertible p^ , since they (and only they) might be transformed reversibly 
by local transformations. The generic orbit of pure states equivalent by local unitary operations has iV 2 — N — 1 
dimensions, and this number decreases in the case of a degeneracy in the vector of the Schmidt coefficients 
On the other hand, the 'noncomparable' states cannot be joined by a local transformation in any direction. The 
existence of such pairs of states (analogous to the space-like events in special relativity) was recently discussed in 
|p6|,[7 U6^j63]| . In fact, Nielsen conjectured that the probability of picking at random two incomparable states out of 
the set of all N x N pure states (according to the natural, rotationally invariant measure) tends to unity in the limit 
N — ► oo. We may now support this conjecture with a simple geometrical argument: the states close to the faces of 
the Schmidt simplex are incomparable, (see Fig. 6c) and it is well known that for large dimensional simplices almost 
all its volume is contained in a thin 'skin' close to its faces. Furthermore, the natural, unitarily invariant measure on 
(£pjv-i g enera t es the Hilbert-Schmidt measure in the Schmidt simplex Jl7 48 1, according to which the degeneracies 
in A are avoided: the probability distribution exhibits maxima in vicinity of the faces of the simplex, while its minima 
occur at the center of the simplex and along the symmetry planes (bisetrices of the triangle) including the corners. 



E. Probabilistic transformations of pure states 



Finally we discuss the case when the majorization condition (3.8) is not fulfilled, so that deterministic local trans- 
formation of the pure state into \(j>) is not possible. It is not possible to create entanglement by means of LOCC if 
the initial state is separable, but it can be shown that if the numbers of nonzero components in both Schmidt vectors 
are the same then one may still transform into \4>) with a non-zero probability P of success. More precisely it was 
shown by Vidal [^6| that the optimal protocol yields 

Pm _ m = g| feA ' W = , nun, f#. (3.10) 



ke[Cm E 4 =fc A i(0) ~ fc efi>] E k ((f>)' 



This statement may be rephrased by saing that the optimal probability P(\ip) — > \<fi}) is the greatest weight p such 
that the submajorization relation holds [p3[ , 

A^^pA^. (3.11) 

Probabilistic transformations are important because they make entanglement purification procedures possible. 

Figure 7 shows the probability of accessing different regions of the Schmidt simplex for pure states of a 3 x 3 system 
for four different initial states The shape of the black figure (p — 1 represents deterministic transformations) is 
identical with the set 'Future' in Fig. 6. The more entangled final state \4>) (closer to the maximally entangled state 
in the center of the triangle), the smaller probability p of a successful transformation. Observe that the contour lines 
(plotted at p = 0.2, 0.4, 0.6 and 0.8) are constructed from the iso-entropy lines H a for a — > and a — > oo (see Fig. 1.) 
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FIG. 7. Probabilistic transformation of pure states in the simplex of Schmidt coefficients for 3x3 pure states: white 
(x) represents the initial state \tjj} (+) the interconvertible states, and black (*) the maximally entangled state. The gray 
scale represents the probability p, with which a given state may be obtained by local transformations from Initial state 
\ifi) is characterized by the Schmidt coefficients A = (Ai,A2,l — Ai — A2), where (Ai,A2) = (a) (0.7,0.25), (b) (0.6,0.27), (c) 
(0.45, 0.45), and (d) (0.8, 0.1). Due to the degeneracy in A in the last two cases there exist only three different interconverible 
states. 



IV. CONCLUDING REMARKS 



We presented a brief review of various measures of entanglement and ways to define them. Analyzing geomet- 
rical properties of the pure states entanglement we have shown that there exist several non-equivalent measures of 
entanglement which generate different ordering of states. Four different measures often discussed in the literature: 
robustness of entanglement (negativity or cross norm), entropy of entanglement, Bures distance to the closest mixed 
state and the distance to the closest pure state are functions of the generalized Renyi-like entropies of the vector of 
the Schmidt coefficients (2.4) with a = 1/2, 1, 2 and a — > 00, respectively, - see Table 1. 
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Measure of 
entanglement 


Formula for 
pure states 


Monotone function 
of H a with 


Schmidt rank 


number of nonzero A; 


a = 


Maximal fidelity 


F max = exp(H 1/2 )/N 


a=\ 


Negativity 


t = exp(# 1/2 ) - 1 


Robustness 


jVF maj£ -l = exp(ff 1/2 )-l = t 

. I : 


Cross-norm 




Entropy of entanglement 


E = Hn = - YliLi m \ 


a = 1 


Relative entropy 


minsBp S(p\\p s ) = E 


Entropy of formation 


E F = E 


Entropy of distillation 


E D = E 


Minimal Bures distance 

to a separable (mixed) state 


rmn SE p{D B ) = 2[1 - exp(-iff 2 )] 


a = 2 


Generalized concurrence 




Minimal distance to 

a separable (pure) state 




OL — ► OO 


Fubini-Study distance 


mm S Ep(D FS ) = arccos[\/A niax j 


trace distance 


rains Ep(Dtr) =2y/l- A max 


Hilbert-Schmidt distance 


mm SEP (D HS ) = y/2(l- A max ) 


Bures distance 


min s ep(Db) = V^i 1 - V^max) 



Table. 1. Measures of the entanglement for N x N bi-partite system ordered according to the properties at the 
manifold of pure states. 



Out of these measures only the entropy of entanglement is additive, if pure states are concerned. The remaining 
quantities (not equivalent to S) are not, but each measure may be related with a certain additive generalized en- 
tropy H a . Thus no single measure seems to be a priori distinguished. On the other hand, some of them may be 
preferred by certain experimentally feasible protocols of local manipulations of the pure states, (e.g. the distillation 
of entanglement), or by the property (E5') of the asymptotic continuity (8]||. 

The above analysis has direct consequences for the more complicated problem of the entanglement of mixed states. 
Consider the definition of entanglement of formation of a mixed state || 

L 

EoF(p) := min^p^(l^)), Pi > 0, (4.1) 

i=l 

where the minimum is taken with over all pure states decompositions p — Y]f—i Pi\ip) (ip\ of a finite length L, normalized 
by ^2f—iPi = 1. In this expression one makes use of the entropy of entanglement, E(\tp)) = Hn(\i())), but one may 
also apply other generalized entropies H a (\ip)) to generate reasonable measures of the mixed states entanglement 
It is easy to see that these generalized entropies of formation induce different orderings in the set of density matrices. 

Time evolution of pure states under such transformations reveals an analogy to special relativity: We first divide 
the set of pure states into equivalence classes consisting of the orbits under local unitary transformations. Then we 
find that for each point in this space of equivalence classes of pure states there exists a 'light cone' which divides the 
space into the 'Future' (F), the 'Past' (P) and the non-comparable states (C). In the Schmidt simplex the role of the 
light cone is played by the accessible region limited by (hyper)faces spanned by any set of (N — 2) edges determined 
by the T-transforms. Thus it is a generalized cone bounded by hyperplanes. This "causal structure" (to pursue the 
analogy to relativity) induces a partial ordering of the states but there is no complete ordering and in this sense then 
there does not exist a preferred measure of entanglement, just as there does not exist an absolute time in relativity. 
On the other hand (unlike in special relativity) the underlying "spacetime" is not maximally symmetric, although it 
comes close to be so. In such situations one often finds that one time slicing is more "natural" than another. Pursuing 
the analogy further one expects that a particular measure of entanglement can be singled out by the physical context. 
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cussions. We appreciate helpful remarks raised by the referee. Financial support by Komitet Badah Naukowych in 
Warsaw under the grant 2P03B-072 19, the NFR, and the European Science Foundation are gratefully acknowledged. 
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APPENDIX A: DISTANCES IN WEYL CHAMBERS AS DISTANCES BETWEEN ORBITS 



Before discussing possible interpretation of the distance in Weyl chamber let us formulate two simple lemmas related 
to the majorization theory. 

Lemma 1. Let x, y and z be three normalized probability vectors with N components in the descending order 
each, and let z < y. Then the scalar products fulfill 

x ■ z < x ■ y. (Al) 

Proof: For simplicity we start with N = 3 and consider the difference x-y~x-z. It is equal to the sum of three terms, 
xi(yi-zi)+x 2 (y2-Z2)+X3(y3-z 3 ) = (xi-x 2 ){yi-zi) + (x2-x 3 )(y 1 +y2-zi-z 2 )+x 3 (yi+y2+y3~z 1 ~Z2-z 3 ), the first 

two are non-negative and the last one vanishes. In the similar way, for an arbitrary N we represent X^feLi x k(yk — 2fe) 

as X)fcli 1 ( a -fc — x k+x) Si=i(y« — z i)- Due to ordering of the components of x and the majorization relation z -< y this 



expression contains non-negative terms only, and this proves the thesis (Al) □. 

Although the next lemma we shall use is closely related with recent results of Nielsen |58| , we formulate it in the 
form most suitable for our purposes and for consistency provide its proof. 

Lemma 2. Let p be a density matrix, A its spectrum and x denotes the diagonal elements of p represented in a 
certain basis. Then x -< A. 

Proof: Any density matrix may be diagonalized, p = UXU^ . Thus the diagonal elements read Xi = pu = UijXjU^. 
Introducing a unistochastic matrix Dy := |tAj| 2 we may then write x = DX, so due to the Horn lemma |36j ] the 
majorization relation x -< A holds □. 

Now we are ready to formulate statements concerning the Weyl chamber which contains the set of diagonal density 
matrices with all entries ordered. 

Proposition 1: Let h and g be two diagonal density matrices of size N with all components in a descending order, 
so they belong to the same Weyl chamber. Let U and V denote arbitary unitary matrices of size N. Then 

a) D HS (h,g) < D HS {UhU\VgV^), 

b) the Euclidean distance between any two points in the Weyl chamber is equal to the Hilbcrt-Schmidt distance 
between the global orbits they generate. 

Proof:. Let W = U'V. Applying the definition [Djjs(h, <?)] 2 = Tr(ft — g) 2 we may rewrite the difference of the 
squared distances, Q = [D HS {UhU\ VgV^)] 2 - [D HS (h,g)} 2 as -2TxW^hWg + 2Txhg. Let h! denote the vector of 
the diagonal elements of W^hW, so the difference we compute is Q — 2(h ■ g — h! ■ g). Applying lemma 2 we obtain 
h! -< d, and using lemma 1 we conclude that the difference Q is non-negative, which gives the statement a) . Thus the 
HS-distance between two orbits, UhU^ and VgV\ equal to the HS-distance between closest points h and g, reads 
Dus{h,g) = \52,i—Ahi — gi)] 1 / 2 . This is just the standard expression for the Euclidean distance between two points 
in the Weyl chamber which represent both diagonal density matrices, h and g □. 

A similar statement may also be formulated for other distances (trace, Bures), but they do not allow for such a 
simple geometric interpretation as does the Hilbert-Schmidt distance used above. Note that the assumption that both 
density matrices belong to the same Weyl chamber is crucial, since a generic orbit visits the simplex of eigenvalues 
N\ times, corresponding to permutations of the components in the spectrum. 

The above proposition clarifies the geometry of Fig. 5a, which shows the asymmetric part of the simplex of 
eigenvalues of the density matrices. To characterize in an analogous way the part of the Schmidt simplex of pure 
states presented in Fig. 6a, we need to change its geometry. The equilateral triangle from Fig. 6b has to be deformed 
into an octant of the sphere. More precisely, we take the vector of the Schmidt coefficients A = {Ai, Aat} of a pure 
state of a N x N bi-partite system and transform it to {\\, A^}. Then the N — 1 dimensional Schmidt simplex is 
mapped into a 1/2 N+1 part of the hypersphere S' JV_1 , equipped with a natural Riemannian distance (for = 3 it is 
just the length of the arc at the sphere S 2 ). This geometry allows us to formulate an analogous 

Proposition 2: Let A and p be two vectors of Schmidt coefficients, representing two pure states \<j>) and \tp) of 
a N x N bi-partite system. Assume that all their components are ordered in a descending manner, so both points 
belong to the same Weyl chamber. Let U = U\®U2 and V = Vx ® V 2 denote arbitary local unitary transformations. 
Then 

a) D FS (\t/j), |0)) < Dps(Ux ® U 2 \il>),Vi ® V 2 \4>)), 

b) the Riemannian distance between any two points in the Weyl chamber (a subset of the deformed Schmidt 
simplex) is equal to the Fubini-Study distance between the local orbits they generate. 

Proof:. In the analogy to the proposition 1, the proof of this proposition may be based on the lemmas 1 and 2. 
Alternatively, one can make use of lemma 1 of a recent paper of Vidal et al. (p3[, which states that the maximal 
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fidelity F — \{(j>\U\ ® l^lV')! 2 is obtained for U\ = U2 = I and is equal to F max = (J2k=i \AfcMfc) 2 - This fact allows 
us to establish part b) of the proposition, since the Ricmannian distance between two points in the Weyl chamber 
Dr(\, fl) is equal to cos _1 (VF max ) = i cos _1 (2F max — 1) and coincides with the Fubini-Study distance between pure 
states \4>) and \ip) □. 

Note that the both propositions resemble the situation for a fibre bundle, where a metric on the total bundle space 
gives a shortest distance between the fibres that can then be used to define a metric on the base manifold. However, in 
our situation not all orbits need to be the same: the dimensionality of local orbits generated by degenerated Schmidt 
vectors (located at the boundary of the Weyl chamber) is analyzed in ref [Q . 
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